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Abstract — Though cooperative relaying is believed to be a 
promising technology to improve the energy efficiency of cellular 
networks, the relays' static power consumption might worsen 
the energy efficiency therefore can not be neglected. In this 
paper, we focus on whether and how the energy efficiency of 
cellular networks can be improved via relays. Based on the 
spatial Poisson point process, an analytical model is proposed to 
evaluate the energy efficiency of relay-assisted cellular networks. 
With the aid of the technical tools of stochastic geometry, we 
derive the distributions of signal-to-interference-plus-noise ratios 
(SINRs) and mean achievable rates of both non-cooperative 
users and cooperative users. The energy efficiency measured by 
"bps/Hz/W" is expressed subsequently. These established expres- 
sions are amenable to numerical evaluation and corroborated by 
simulation results. 

I. Introduction 

With the rapid proliferation of user equipments (UEs) and 
base stations (BSs), the energy consumption of cellular net- 
works grows quite amazing. It is shown in [ 1] that the energy 
consumed by BSs accounts for nearly 60-70% of the total 
network energy. Therefore, improving the energy efficiency 
of BSs becomes the fundamental challenge to realize green 
communications. 

Recently, cooperative relaying has been in-depth studied as 
a promising way to reduce the energy consumption of cellular 
networks 12[]|l3|], which can extend coverage, increase capac- 
ity and provide flexible cost-effective deployment options as 
well 01 . Compared with BSs, RSs cover a much smaller 
area and require lower transmit power. Meanwhile, the UEs 
covered by RSs mostly enjoy much higher average signal-to- 
interference-plus-noise ratios (SINRs). As RSs do not have 
a wired backhaul connection, their deployment cost is largely 
reduced. Consequently, using cooperative relaying technique is 
capable of lowering the transmit power consumption in cellular 
networks and the low power RSs are easy to be deployed 
without modifying the current cellular infrastructure. 

Generally, the total power consumption of both BSs and 
RSs consists of two parts, i.e., the transmit power consump- 
tion and the static power consumption. Compared with the 
transmit power consumption, the static power consumption 
which is contributed by signal processing, battery backup, 
site cooling and etc., occupies a substantial part of the whole 
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power consumption. Therefore, the transmit power saved by 
cooperative relaying technique may not compensate for the 
additional power consumption of RSs which leads to higher 
power consumption in cellular networks. 

The aforementioned discussions motivate us to investigate 
whether and how the energy efficiency of cellular networks 
can be improved via relays. The related works mostly focus 
on regular network deployments, such as the hexagonal grid 
model, which are merely estimated by Monte Carlo simu- 
lations [2] [5] [6]. However, the spatial deployment observed 
in actual communication networks is usually far from being 
regular. In addition, the interference and the signal strength 
at a receiver critically depend on the spatial positions of the 
interfering transmitters, thereby mathematical techniques are 
needed to explicitly model the node distribution. Stochastic 
geometry theory has recently emerged as an essential tool to 
model and quantify interference, mean achievable rate, and 
coverage in cellular networks which are verified to be close 
to the actual networks [7]. Consequently, this paper proposes 
an analytical model to evaluate the energy efficiency of relay- 
assisted cellular networks using stochastic geometry, instead 
of relying on system simulations only. The distributions of 
SINRs and mean achievable rates of both non-cooperative 
UEs and cooperative UEs are derived. We use the ratio of 
spectral efficiency to power consumption with unit bps/Hz/W 
as the metric of energy efficiency and further present the 
theoretical expression. Simulation results present how the 
intensity of users affects the energy efficiency and how much 
RS's static power consumption should be controlled below 
in order to obtain some energy saving gains. Notably, our 
work serves as a pioneering effort on the network design and 
planning, especially with respect to energy efficient wireless 
communications . 

II. Network Model and Metrics 

We consider a relay-assisted downlink cellular network 
shown in Fig. Q] The locations of BSs are modeled ac- 
cording to a homogeneous Poisson point process (PPP) $t 
with intensity Ah, and RSs are located according to another 
homogeneous PPP $ r with intensity A r . Specifically, for a 
given PPP, the number of points in a bounded area is a Poisson- 
distributed random variable and the points are uniformly- 
distributed within the area. A realization of the relay-assisted 
cellular network with Poisson distributed BSs and relays 
portrayed by a Voronoi tessellation is given in Fig. [2] Here, 




Fig. 1. Downlink relay-assisted cellular network model. 
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Fig. 2. A Voronoi network topology, where each Voronoi cell is the coverage 
area of a BS and each small circle represents a relay coverage area. 

we use the term cell to refer to a Voronoi region with a 
random area |@] of each BS, which, as noted in [9], more 
closely reflects actual deployments that are highly non-regular. 
We denote the area of the Voronoi cell by S. There is 
no known closed-form expression of the probability density 
function (PDF) of S, whereas a simple approximation lllOll 
has been proved sufficiently accurate for practical purposes. 
Considering scaling, the approximate PDF of the size of a 
single-cell coverage area is given by 

= TT\/?( SA &)* exp(-^A 6 )A 6 (1) 
15 V 2/k I 

Each RS has a circular coverage with radius R r . The coopera- 
tive UEs are distributed according to a homogeneous PPP with 
intensity A c in each circular area, and the non-cooperative UEs 
are distributed according to another homogeneous PPP with 
intensity A„ c in the whole network. All the PPPs are mutually 
independent. 

We assume that each UE is served by its closest base station 
which is at a distance r\g. Since the BSs are distributed as 
PPPs, it follows that is Rayleigh distributed with parameter 
l/^fhFXb [11]. If the UE is also located in a RS's coverage, 
the distance tr between the RS and the UE follows the 
distribution = 2r/R 2 .. The transmit power is Pb for all 

the BSs and Pr for all the RSs. For the sake of convenience, 
we adopt a standard path loss propagation model with path 



loss exponent a > 2. Regarding fading, we assume that the 
link between the serving access point (either a BS or a RS) and 
the served UE experiences Rayleigh fading with parameter fi. 
The set of interfering entities is the ones that use the same 
subchannel as the served UE and the fading of interference 
links may follow an arbitrary probability distribution, which 
is denoted by g. For all receivers, the noise power is a 2 . 

The different PPPs corresponding to different entities in 
the network interact in a complicated way, which makes a 
rigorous statistical analysis extremely difficult. For example, 
a cooperative UE may be covered by more than one RS that 
leads to the delicate issue of relay selection, and furthermore 
renders the distributions of subchannel usage among RSs 
and BSs intrinsically correlated. To overcome the technical 
difficulties due to spatial interactions, we propose a two- 
scale approximation of the network model for the subsequent 
analysis motivated by the fact that the covered area of a 
RS is significantly smaller than that of a BS. The two-scale 
approximation consists of two views, i.e., the macro-scale view 
and the micro-scale view. The macro-scale view concerns an 
observer outside the coverage area of a RS, and the whole 
coverage area of the RS shrinks to a single point, marked 
by the number of cooperative UEs. The micro-scale view 
concerns an observer inside the coverage area of a RS, and 
the coverage area is still a circular with radius R r where the 
cooperative UEs are spatially distributed. By such a two-scale 
approximation, a cooperative UE can only be covered by a 
unique RS, and the coverage area of a RS can only be within 
a unique Voronoi cell of a BS. 

A. Distribution of UEs 

Let U c be the number of cooperative UEs in a RS's cov- 
erage and from our model we have U c ~ Poisson(X c TrR 2 ). 
The number of non-cooperative users served only by BSs is 
denoted by U nc , and it is characterized as follows. Let S be 
an arbitrary single-cell coverage area, then the number of non- 
cooperative users in this cell is a Poisson random variable with 
mean X nc S. Conditioned upon S, the probability generating 
function (PGF) of the unconditioned U nc is thus given by 

/■OO 

G(z) = I exp(A nc (z - l)S)f(S)dS 

where f(S) is given by (Q]i, and the distribution of U nc is 
therefore given by the derivatives of G{z), 

P{U nc = i} = — ±-L, » = 0,1,... (3) 

Similarly, we can obtain P{U r — i}, where U r represents the 
number of relays in a Voronoi cell with mean \ r S. 

B. Spectrum Allocation and Usage 

The available spectrum for a cell is evenly divided into 
M subchannels. We assume that RSs are given a separate 
frequency band from BSs, such that there is no cross-tier 



interference. Thus a partition (Mb,M r ) is made, where Mb 
subchannels are available for BS transmissions and M r = 
M—Mb subchannels are available for RS transmissions. When 
the current UEs are more than the number of subchannels 
in a serving access point, they are served by time-sharing 
with equal time proportion. To distinguish between non- 
cooperative UEs and cooperative UEs for BS, Mb are further 
divided into two parts, i.e., Mb 1 and M(, 2 , which are used for 
non-cooperative transmissions and cooperative transmissions, 
respectively. From the macro-scale view, for a BS, the access 
required from the cooperative UEs is equivalent to that from 
the RS. In this case, we divide Mb 2 into p parts, where 
p = Mb 2 /M r , and due to the broadcast nature of wireless 
medium, each part of the subchannels are used for both BS- 
to-RS link and BS-to-UE links. 

The subchannels in each partition {Mb 1 , Mb 2 , or M r ) are 
uniformly and independently selected by BS or RSs, and it 
suffices to analyze an arbitrary one of them. Firstly, let us 
examine the probability that a given subchannel in Mb 1 is 
used by a BS, which can be obtained via normalizing the 
average number of subchannels used by the total number of 
subchannels Mb 1 - The same goes for Mb 2 and M r given as 
follows. 

PbusyM =j^j-J2 min ^' M bl )P{U nc = i} (4) 
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min(i, p)P{U r = i} 
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busy.r = Tf~yi min (^ M r )P{U c = i} (6) 
r i=0 

The spatial process of RSs that use a given subchannel is 
the independent thinning of the original PPP of RSs $ r by 
the probability Pbusy,r, denoted by $ r with intensity A,, = 
X r Pbusy,r- The term "independent thinning" means that $ r 
can be viewed as obtained from <i> r by independently removing 
points with probability l—Pbusy,r- Similarly, the spatial process 
of BSs that use a given subchannel (either in Mb 1 or Mb 2 ) is 
the independent thinning of the original PPP of BSs by the 
probability PbusyM ( or Pbusy,b 2 ), denoted by <£> b with intensity 
\' b = XbPbusyM (or X' b = X b Pbusy,b 2 )- These two independently 
thinned PPPs will prove useful in the subsequent analysis. 

C. Power Model 

We model the relation between the average power consump- 
tion and the average radiated power per site in a linear fashion 

P i =a i P i + b i , fori = B,R (7) 

where Pi and Pi denote the average consumed power and the 
radiated power per site, respectively, while B, R represent 
BS and relay. The coefficients accounts for the power 
consumption that scales with the average radiated power due 
to amplifier and feeder losses as well as cooling of sites. 
The terms bi denotes power offsets which are consumed 
independently of the average transmit power. These offsets 



are, amongst others, due to signal processing, battery backup, 
as well as site cooling jfj. 

D. Energy-Efficiency Function 

The energy-efficiency evaluation function of cellular net- 
works in this paper is defined as the ratio of spectral efficiency 
to power consumption with unit bps/Hz/W. Here we use 
the mean achievable rate to measure the spectral efficiency. 
According to our model, the cellular network is homogeneous 
and the energy efficiency of the whole network is equivalent 
to the energy efficiency of an arbitrary Voronoi cell which can 
be expressed as 



Q 



(8) 



Pb + NP R 

where r denotes the cell's total mean achievable rate provided 
by both BS with consumed power Pg and RSs with consumed 
power NPr, and N is the average number of RSs in the cell. 

III. Analysis of Mean Achievable Rate and Energy 
Efficiency 

In this section, we derive the mean achievable rates of 
cooperative UEs and non-cooperative UEs at first, and then 
give the energy efficiency expression based on the definition 
in lll-DI For each type of UEs, we begin with general settings, 
and then simplify the general results under several specific 
parameters. The mean achievable rates are the averaged instan- 
taneous achievable rates over both channel fading and spatial 
distributions of UEs and access points. 

A. Non-Cooperative UEs 

1) General Case: The UE is served by occupying one of 
subchannels with transmit power Pb = jf- in the BS at a 
distance r, g. is the fading of an interference link, and R. is 
the distance between the UE and an interfering access point, 
thus the SINR experienced by the UE is given by 7 = Pbhr 



where exp(p), Ib — ^ ie ^'^ bQ PbgiRi a is the interference 
from the "BS" tier (excluding the serving BS itself which is 
denoted by bo). We can thus derive the CDF of 7 as 



Pr( 7 < T) 
1 - / Pr 

1 - 
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Pr( 7 > T) 
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(9) 



where f(r) = 2n\bre x >> 7vr " . By the definition of Laplace 
transform we get 
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exp(-2TT\ b (1 - C g (sv- a P b ))vdv) (10) 



where the last equality follows from the PGF of PPP 111 ill . By 
exchanging the integral order, we obtain 



(A) 



r ^J^ Eg \ g ifr(-l, sPbg r--)-T(-^))} (11) 
2 a L V a a /J 

Then we can obtain the CDF of 7 by substituting ( TTOb and 
CO) into ©. 

Now we evaluate the mean achievable rate. Since for a 
positive random variable X, E[X] = J t>Q P{X > t}dt, we 
have 



T nc = E[ln(l + 7)] 

/„ ./„ V /'.'• " / 

^* 2 ( c t-i) /nr a (e t -l)\ 

e —C h ^ % 1) )dtf(r)dr 
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dtdv 



(12) 



where <^(T, a) is given by 



<p{T,a) = -(nT)-E g 



g i( T (—,»Tg)-r(—)) (13) 

The integrals in (fT2l can be evaluated by numerical methods, 
while they can be simplified to a concise form in the special 
case. 

2) Special Case: When the interference experiences 
Rayleigh fading with mean /1, i.e., g ~ exp(p), and a = 4. 
The results are simplified as follows. 
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2 / ucr 2 (e t - 1) 
y2 / 2 dy, and k(T) is given by 



dt (15) 



k(T) = ir\ b (l +Pi rasJ/ ,6 1 V / T^ - arctan-^=J ) (16) 



S. Cooperative UEs 

1 ) General Case: To obtain the mean achievable rate of a 
cooperative UE, we have to derive the distributions of SINRs 
(7o> 7ii 72) and mean achievable rates (to, t\, T2) for BS-to- 
UE link, BS-to-RS link and RS-to-UE link, respectively. Based 
on our two-scale approximation, we use 71 to approximate 70 
in the subsequent analysis. Fig. [3] in [IV] verifies the reasonable- 
ness of this approximation. The transmission scheme termed 



coded cooperation with selective source-to-destination parity 
transmission lfl2ll is introduced here for cooperative UEs. The 
scheme has two modes depending on whether the RS decodes 
the message from BS correctly. Accordingly, different modes 
have different mean achievable rates. The specific analysis is 
presented below. 

The probability of correctly decoding for the RS is denoted 
by Pd and expressed as 



P d = Pr(7i > T th ) = T:\ b J exp ( - 
+nX b v( 1 - ip(T th ,a) 



pvia 2 T th 



P, 



1 



dv (17) 



Pbusy ,62 ' 

where T t / l =2 i?th /' 3 — 1, R t h is the target rate of RS correctly 
decoding and /3 is the cooperative parameter in coded coop- 
eration. 

a) Mode 1: When the RS fails to decode, the cooperation 
is degraded to a direct transmission. 

{1 T > T 

W( 71 <T) T ~ T th (18) 
iM-n<T«0' 1 < - ±th 

The conditional PDF of 71 denoted by / 7l (i | 71 <Tth) is 
obtained by differentiating ( TT~8T > with respect to T. Then we 
get the mean achievable rate in this mode as 



ln(l + t)/ 71 (t 1 71 <Ttfc)dt 



ln(l+T th )- 



1 



' tfe Pr(7i < t) 



l + t 



dt (19) 



b) Mode 2: When the RS decodes correctly, the data is 
transmitted from BS to UE and the parity is transmitted from 
RS to UE. 

For BS-to-RS link, the derivation of n is similar to non- 
cooperative UEs in IIII-AI by just revising Pbusy. b x to Pbusy, b 2 - 
For BS-to-UE link, 

f 0, T < T th 

Pr{ii<T\"f X >T th ) = I pr( 7l <r)-Pr( 7l <r fh ) T>T 

I l-Pr( 7 i<T th ) ' — th 

(20) 

The conditional PDF of jx denoted by f 11 (f | 71 > T tb ) is 
obtained by differentiating d20b with respect to T. Then we 
get the mean achievable rate of BS-to-UE link as follows 

/•OO 

r = / ln(l + t)f yi (t\ji>T th )dt 

J T t h 

= -^(n - (1 - P d )r ml ) (21) 

For RS-to-UE link, the UE is served by occupying one of 
subchannels with transmit power P r = jf* in the RS at a 
distance r. Denote by I r the interference strengths from the 
"RS" tier. The CDF of 72 is given by 

Pr{ l2 < T) = 1 - e-^^C^i^^lf^dr (22) 

JO "r 



Following the same procedure as in the derivation of (10), we 
get the Laplace transform of as 

Ci r (s) = exp (-ir\' r (sP r )*T(l - ^)E a (<^)J (23) 

Substituting (l23l into d22l with v = r 2 , we get the distribution 
of 72. Similar to ( TTZl ). r 2 is given by 

(24) 

where VK**) i s given by 

V>(a) = -7r^A' r r(--)E 9 ( 5 *) (25) 
a a 

According to the capacity formula in | fl2t] . we get the mean 
achievable rate in this mode below 

Tm2 = /3t + (1 - P)t 2 (26) 

Integrating r m i and r m 2, the mean achievable rate for a 
cooperative UE is obtained as 

r c = (1 - P d )r m i + iV m2 (27) 



2) Special Case: When the interference experiences 
Rayleigh fading and a = 4, the results are simplified as 




C. Energy Efficiency 



The energy efficiency in our network is expressed as 

Q= I sl+T t (31) 
Pb + NP R 

where r sl and t s2 are given by 

+oo 

Tsl = ^ mia ( M bl ) i)?{Unc = i}r„c (32) 
i=l 

+OO+OO 

r S 2 = y^y^ min(p, i)min(M r , j)P{U r = i}P{U c = j}r c 

i=lj=l 

(33) 



IV. Numerical Results 

In this section, we present comparison of the Monte Carlo 
simulations and the analytical results to evaluate the energy 
efficiency of a downlink relay-assisted cellular network. The 
default configurations of system model are shown in Table Q] 

Figure |3]displays the SINR distributions of cooperative UEs 
and non-cooperative UEs. The curve of the non-cooperative 
UE reveals that the simulation result matches the analytical 
result well, thus corroborating the accuracy of our theoretical 
analysis. For the cooperative UE, the theoretical curve denotes 
the SINR distribution of BS-to-RS link, which verifies that the 
SINR distribution of BS-to-UE link can be approximated by 
that of BS-to-RS link. 

Figure|4]illustrates how the energy efficiency of this network 
changes with various X nc and A c . We observe that the energy 
efficiency Q increases with A nc and A c until reaching Q op t with 
value 0.78bps/Hz/W and then tends to keep constant after a 
slight decrease. The reason should be that in the early stage 
as the intensities of both cooperative UEs and non-cooperative 
UEs begin to increase, the idle subchannels are sufficient for 
the existing UEs which results in the increase of the energy 
efficiency. After reaching Q op t, more subchannels become 
occupied by BSs and RSs which incurs more "BS" and "RS" 
tier interference and leads to the decrease of energy efficiency. 
But soon after the intensity of UEs increases sufficiently large, 
all the subchannels are persistently occupied by BSs and RSs 
with the UEs served by time-sharing, and then the interference 
saturates, which makes the energy efficiency unchanging. 

Figure [5] shows how the offset power bn and the Tx power 
of RS affect the energy efficiency. We use the variable r] to 
change the value of bn beginning with the default value 32. 0W, 
i.e., bji = 7/x32.0W. The curves reveal that the simulation 
results match the analytical results well which corroborate 
the accuracy of our theoretical analysis of energy efficiency. 
The figure shows that the energy efficiency of the RS with 
24dBm Tx power case outperforms that of the 33dBm and 
38dBm cases, and the energy efficiency decreases with the 
increasing of rj. In addition, we get the significant finding here 
that there exist thresholds of RS's offset power to compare 
the energy efficiency metric in the cellular networks with 
and without relays, which verifies the research motivation of 



TABLE I 
System Parameters 



Symbol 


Description 


Value 


A 6 


density of BSs 


10-" BS/m 2 




density of RSs 


9 x ltr : ' RS/m^ 


Pb/Pr 


transmit power 


43dBm/33dBm (BS/RS) 


M 


number of subchannels 


300 (M b =285, A/ r =15) 


Rr 


radius of RS 


20m 


a 


pathloss exponent 


4 




Rayleigh fading parameters 


1 (normalized) 


a 2 


noise power 


-80dBm 


P 


cooperative parameter 


0.6 


Rth 


target rate of RS correctly decoding 


0.5bit/s/Hz 


a B ,b B 


power model parameters of BSs 


22.6, 412.4W \6] 


a-R, b R 


power model parameters of RSs 


5.5, 32.0W \6] 




Fig. 3. CDFs of SINRs for non-cooperative UEs and cooperative UEs. 

Energy Efficeincy vs. Intensity ot UE 




Fig. 4. Energy efficiency versus the intensities of non-cooperative UEs and 
cooperative UEs. 

this paper. Specifically, if the offset power of RS is smaller 
than the threshold, then it is worthwhile to introduce relays 
into the cellular networks, otherwise not. For example, when 
the Tx power of RS is 24dBm, the thresholds in this case 
is r]th = 0.95, i.e., = 30AW, then the relay-assisted 
networks will have better energy efficiency when b R < b 1 ^ 
compared with conventional networks without RSs. 

V. Conclusions 

In this paper, we proposed an analytical model for evaluat- 
ing the energy efficiency of a relay-assisted cellular network 
using stochastic geometry approach instead of merely relying 
on system simulations. We introduced coded cooperation as 
the transmission scheme, and derived the distribution of SINRs 
and mean achievable rates of both non-cooperative UEs and 
cooperative UEs. And then the analytical expression of energy 
efficiency for our network is obtained. In our simulations, a 
3D surface of energy efficiency was presented to reveal the 
relationship among the energy efficiency and the intensity 
of each type of UEs. Besides, we observed that the static 
power consumption of RSs behaves as an important factor 
to determine whether cooperative relaying can save energy of 
cellular networks or not. For certain network configurations, 



Fig. 5. Energy efficiency versus the offset power bft and the Tx power of 
RS. 

there exist static power thresholds of RSs, under which the 
energy efficiency of the cellular networks can be meliorative. 
The work in this paper sheds a useful insight into the real 
cellular network design and planning. 
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